Solar dynamo as host power pacemaker of the Earth global climate 
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It is known that the so-called problem of solar power pacemaker related to possible existence 
of some hidden but key mechanism of energy influence of the Sun on fundamental geophysical 
processes is one of the principal and puzzling problems of modern climatology. The "tracks" of 
this mechanism have been shown up in different problems of solar-terrestrial physics for a long time 
and, in particular, in climatology, where the solar-climate variability is stably observed. However, 
the mechanisms by which small changes in the Sun's energy (solar irradiance or insolation) output 
during the solar cycle can cause change in the weather and climate are still unknown. 

We analyze possible causes of the solar-climate variability concentrating one's attention on the 
physical substantiation of strong correlation between the temporal variations of magnetic flux of 
the solar tachocline zone and the Earth magnetic field (Y-component). We propose an effective 
mechanism of solar dynamo-geodynamo connection which plays the role of the solar power pacemaker 
of the Earth global climate. 

Keywords: Solar dynamo-gcodynamo connection; Earth's global climate; Solar- climate variability; Solar 



I. INTRODUCTION 

It is known that in spite of a long history the nature of 
the energy source maintaining a convection in the liquid 
core of the Earth or, more exactly, the mechanism of the 
magnctohydrodynamic dynamo (MHD) generating the 
magnetic field of the Earth still has no clear and unam- 
biguous physical interpretation (see [1] and refs. therein). 
The problem is aggravated because of the fact that none 
of candidates for an energy source of the Earth magnetic- 
field [1] (secular cooling due to the heat transfer from the 
core to the mantle, internal heating by radiogenic iso- 
topes, e.g., 40K, latent heat due to the inner core solid- 
ification, compositional buoyancy due to the ejection of 
light clement at the inner core surface) can't in principle 
explain one of the most remarkable and mystic phenom- 
ena in solar-terrestrial physics, which consists in strong 
(negative) correlation [2] [3] between temporal variations 
of the magnetic flux in the tachocline zone (the bottom 
of the Sun convective zone) [4] and the Earth magnetic 
field [5] (Figure 1). 

At the same time, supposing that the transversal (ra- 
dial) surface area of tachocline zone, through which a 
magnetic flux passes, is constant in the first approxima- 
tion, we can consider that the magnetic flux variations 
describe also the magnetic field temporal variations in the 
tachocline zone of the Sun. In this sense, it is obvious 
that a future candidate for an energy source of the Earth 
magnetic field must play not only the role of a natural 
trigger of solar-terrestrial connection, but also directly 
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generate the solar-terrestrial magnetic correlation by its 
own participation. 

The fact that the solar-terrestrial magnetic correlation 
has, undoubtedly, fundamental importance for evolution 
of all the Earth's geospheres is confirmed by existence 
of a stable and strong correlation between the temporal 
variations of the Earth magnetic field, the Earth angular 
velocity, the average global ocean level and the number of 
large earthquakes (with the magnitude M>7), whose gen- 
eration is apparently predetermined by a common phys- 
ical cause of unknown nature (see Figure 1). 

On the other hand, it is clear that understanding of 
the mechanism of solar-terrestrial magnetic correlation 
can become the clue of so-called problem of solar power 
pacemaker related to possible existence of some hidden 
but key mechanism of energy influence of the Sun on the 
fundamental geophysical processes. It is interesting, that 
the " tracks" of this mechanism have been observed for a 
long time and manifest themselves in different problems 
of solar-terrestrial physics and, in particular, in clima- 
tology, where the mechanisms by which small changes in 
the Sun's energy output during the solar cycle can cause 
change in the weather and climate have been a puzzle 
and the subject of intense research in recent decades. 
Thus it becomes obvious that purposeful or unpurposed 
neglect of the mechanism of solar power pacemaker in 
any point or multi-zonal models of the Earth global cli- 
mate can result in serious errors in interpretation of the 
experimental temperature and other geophysical trends, 
especially in the compilation of short-term and, all the 
more, long-term forecasts. 

In this paper we consider hypothetical particles ( 57 Fe 
solar axions) as the main carriers of the solar-terrestrial 
connection, which can transform into photons in exter- 
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nal fluctuating electric or magnetic fields by virtue of 
the inverse coherent Primakoff effect [6]. At the same 
time we ground and develop the axion mechanism of so- 
lar dynamo - geodynamo connection, where the energy 
of axions is modulated at first by the magnetic field of 
the solar tachocline zone (due to the inverse coherent Pri- 
makoff effect), and after that is resonantly absorbed in 
the iron core of the Earth, thereby playing the role of an 
energy source and modulator of the Earth magnetic field. 
Justification of the axion mechanism of solar dynamo - 
geodynamo connection and its account within the frame- 
work of the bifurcation model of the Earth global climate 
on different time scales is the goal of this article. 



consists in determination of the global temperature T(t) 
and its increment AT(t), is reduced to finding the stable 
solution set of equations: 



JLu*{T,t) = T? + a{t)-T t + b{t) = Q, 



where 
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II. PECULIARITIES OF THE BIFURCATION 
MODEL OF THE EARTH GLOBAL CLIMATE ON 
DIFFERENT TIME SCALES 
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As is shown in our papers [10, 11], the basic equation of 
energy-balance model of the Earth global climate is the 
bifurcation equation (with respect to the Earth surface 
temperature (see Figure 2)) of assembly- type catastrophe 
with two governing parameters, which describe insolation 
variations and the Earth magnetic field variations (or the 
variations of cosmic ray intensity in the atmosphere). A 
general bifurcation problem of this energy-balance model 
(see equations (20)-(23) and (26)-(28) in [10, 11]), which 
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U* (T, t) describes with an accuracy up to constant the 
so-called " inertial" power of heat variations in the Earth 
climate system; AU*(T,t) is the variation of U*(T,t); 

is the relative intensity of terrestrial magnetism, T t 
is the average global temperature of the Earth surface 
at the time t, K; AT t is the variation of T t ; So = 
1366.2 Wm~ 2 is "solar constant"; 6 = 0,95 is coeffi- 
cient of gray chromaticity of the Earth surface radia- 
tion; er = 5, 67 • 10 -8 is the Stephen-Boltzmann con- 
stant, Wm~ 2 K- 4 ; r] a = 0, 0295-ftT -1 ; f3 is the accu- 
mulation rate of carbon dioxide in the atmosphere nor- 
malized by unit of temperature, kgK , a M and b^ 
are constants, whose dimensions are Wm~ 2 K~ 2 and 
Wmr 2 K~ Y , respetively; AW(t) is the insolation reduced 
normalized variation; a s is the root-mean-square devia- 
tion; W rec iuced = 5*o + AW(t)<7 s is the reduced annual 
insolation. 

Within the framework of proposed bifurcation model 
(i) comparison of the solution of energy-balance model 
of the Earth global climate and the EPICA Dome C 



and Vostok experimental data of the Earth surface 
palacotcmperature evolution over past 420 and 740 kyr 
is given; (ii) possible sharp warmings of the Dansgaard- 
Oeschger type during the last glacial period due to 
stochastic resonance is theoretically argued; (iii) the con- 
cept of climatic sensitivity of water in the atmosphere, 
whose temperature instability has the form of so-called 
hysteresis loop, is proposed, and based on this concept 
the time series of total fresh water mass (or vice versa 
the global ice volume) over the past 1000 kyr, which 
is in good agreement with the time series of <5 18 con- 
centration in sea sediments (Figure 3), is obtained; (iiii) 
groundlessness of the so-called " CO2 doubling" problem 
is discussed. 

One of the main features of bifurcation model of the 
Earth global climate lies in the wonderful fact that priori 
knowledge of only two governing parameters, which are 
set by the known time series of insolation variations and 
variations of the Earth magnetic field (or the variations 
of cosmic ray intensity in the atmosphere) is required to 
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FIG. 1. Time evolution (a) the variations of magnetic flux in the bottom (tachocline zone) of the Sun convective zone (see Figure 
7f [4]), (b) of the geomagnetic field secular variations (Y-component, nT/year), whose values are obtained at the Eskdalemuir 
observatory (England) [5], where the variations (5Y/St) are directly proportional to the westward drift of magnetic features, (c) 
the variation of the Earth's rotation velocity [7], (d) the variations of the average global ocean level (PDO+AMO, cm/year) [8] 
and (e) the number of large earthquakes (with the magnitude M > 7) [9]. All curves are smoothed by sliding intervals in 5 and 
11 years. The pink area is the prediction region. Note: formation of the second peaks on curves (c)-(e) is mainly predetermined 
by nuclear tests in 1945-1990. 




FIG. 2. The energy fluxes balance on the Earth surface. Here 
G w , G v , Geo are the heat energy power re-emitted 2 by 
liquid water, water vapour and carbon dioxide, respectively, 
G friction is the heat or dissipation energy generated by the 
Earth surface-to-atmosphere bottom layer friction. 



solve the basic equations (l)-(6) of energy-balance model 
(or to determine theoretical temperature trends). Other 
not less interesting feature of this model is the so-called 
principle of structural invariance, which means that the 
shape of global climatic potential (assembly-type catas- 
trophe) 



U(T,t) = -T i + -a(t)T 2 + b(t), (7) 

is structurally invariant on different time scales. In 
other words, the principle of structural invariance of the 
balance equations of climate models evolving on the dif- 
ferent time scales is not only a direct indicator of the 
correctly guessed physics of non-uniformly scaled pro- 
cesses, but it simultaneously specifies unambiguous rules 
for transition from one time scale to the other within the 
framework of one global model as well as for transition 
from the (one-zonal!) model of the Earth global climate 
on any time scale to the multizonal model of global cli- 
mate or weather on a short time scale. It means that the 
system of equations of the multizonal model of global cli- 
mate or weather convoluted into the balance equation of 
one-zonal model must fully keep the structure and prop- 
erties (governing parameters) of the bifurcation model of 
global climate on different time scales. Since the bifurca- 
tion model describes the climatic trends of palcotempcr- 
ature and global ice volume well without consideration 
of the mechanism of solar power pacemaker, the natural 
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FIG. 3. Comparison of the theoretical time series of climatic 
sensitivity A„+„ calculated by equations (3) and (19) from [11] 
(b) with the time series of 8 ls O isotopic concentration (the 
conditional analogue of ice volume) measured in the deep- 
water experiments: (a) Bassinot et al. [12] (solid blue line) 
and Imbrie et al. [13] (dashed red line), (c) Tidemann et al. 
[14] 



question arises here: "Is it possible that this fact contra- 
dicts the assigned task?" Below we will show that there 
is no contradiction here, because actually the mechanism 
of solar power pacemaker is implicitly taken into account 
in the climatic potential (7), and it is non-trivial confir- 
mation of significance of the principle of structural in- 
variance. 



III. AXION MECHANISM OF SOLAR 
DYNAMO - GEODYNAMO CONNECTION 



We have shown that strong correlation between the 
temporal variations of magnetic field of the Earth (Y- 
componcnt) and the magnetic field toroidal component 
of tachocline zone of the Sun really takes place. There- 
upon we have asked ourselves: "May hypothetical solar 
axions[15], which can transform into photons in external 
electric or magnetic field (the so-called inverse PrimakofF 
effect), be the instrument by which the magnetic field of 
the solar tachocline zone modulates the magnetic field of 
the Earth? In other words, may solar axions be an effec- 
tive energy source and modulator of the Earth magnetic 
field?" 

It turns out that it is really possible [6] . Following [6] , 
let us consider without loss of generality the simplified 
axion mechanism of solar dynamo-geodynamo connec- 
tion. As is known, the reaction of the solar cycle that 



produces solar energy is one of axion sources. Since ax- 
ions are pseudoscalar particles they can be emitted in 
nuclear magnetic transitions. On the other hand, since 
the temperature in the center of the Sun is 1.3 keV, some 
nuclei having low-lying nuclear level can be excited ther- 
mally. At the same time, monochromatic axions can be 
emitted in the nuclear magnetic transitions from the first 
thermally excited level to the ground state. Below we 
consider only 14.4 keV solar axions emitted by the Ml 
transition in 57 Fe nuclei, because just these axions can 
be resonantly absorbed in the iron core of the Earth gen- 
erating 14.4 keV 7-quanta by the discharge of the excited 
nuclear level (Figure 4). 

It is interesting, that exactly these 7-quanta with the 
energy 14.4 keV are the supplementary energy source in 
the Earth core, which can pretend to the role of energy 
source of generation and modulator of the Earth mag- 
netic field. At the same time there is a natural question, 
is this energy sufficient for generation of the magnetic 
field of the Earth and how this source can execute the 
role of the modulator of the Earth magnetic field. To 
answer these questions let us briefly consider the axion 
"course of life" inside the Sun before it leaves the Sun. 
It appears [6] that passing through the solar tachocline 
zone (the bottom of the Sun convective zone), where the 
Sun magnetic field is generated, axions can be converted 
into 7-quanta and thereby to decrease the solar axions 
flux to the Earth. As is shown in [6], in this case the 
probability that an axion converts back to a " observable" 
photon inside the magnetic field can be represented by 
the following simple form 



p ~ 

1 07 — 



9ayBL 



(8) 



where g ai ~ 1.64 • 10 9 GeV' 1 is the strength of an axion 
coupling to a photon, L ~ 3.5 • 10~ 7 m is the thickness of 
solar tachocline zone, B ~ 35 T is the conservative value 
for the magnetic field of the active Sun. From this it 
follows that the solar axion flux outgoing beyond the Sun 
is modulated by the value of the Sun magnetic field (see 
(7)). At the same time, it is obvious, that the axion flux 
to the Earth is low during the active Sun and conversely 
it is practically maximal during the quiet Sun. 

Now let us show that the total energy of axions during 
the quiet Sun is sufficient to generate the Earth magnetic 
field. It is not difficult to show [6] that the axion resonant 
absorption rate in the Earth core, which contains the Np 7 e 
nuclei of 57 Fe isotope, is about 



R a ~ 5.2 • 10 
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FIG. 4. Schematic picture of the solar tachocline zone, Earth's liquid outer (red region) and inner (brown region) core. Solar 
axions are resonance absorbed in iron of the Earth core conversing into 7-quanta, which are the supplementary energy source 
in the Earth core. Blue lines on the Sun designate the magnetic field. Note: In the conventional concept, the molten iron of 
liquid phase of Earth's core circulates along a spiraling in columns aligned in the north-south direction, generating electrical 
currents that set up the dipolar magnetic field. The concentration of field lines into anticyclonic vortices (rotating in the same 
as air around a region of high pressure) has been thought to explain the intense magnetic lobes found in Earth's field at the 
top of the core. 



It is known, that the number of Fe nuclei in the 
Earth core is N F 7 e ~ 3 • 10 47 [6] and the average energy 
of 57 Fe solar axions is (E a }= 14.4 keV. If in (9) for an 
axion-nucleon coupling g e JJ ~ 10 -5 [6] to take into ac- 
count the factor 2 related to uncertainty of iron concen- 
tration profile at the Sun, then with an allowance of (9) 
the maximum energy release rate W 7 in the Earth core 
is equal to 

W 1 = R a ■ (E a ) ~ 1 TW. (11) 

Analysis of modern model parameters of thermal state 
of the Earth's core [1] shows that in spite of known dif- 
ficulties in interpretation of the results of evolutionary 
geodynamo simulation, such a thermal power (1 TW) is 
sufficient for generation and maintenance of the Earth 
magnetic field [1]. It is easy to show that it is exactly 
so by the known dependence of magnetic field Be on the 
total ohmic dissipation D in the Earth core 

D ~ 7^ B ^ ( 12 ) 

where 77 is the magnetic diffusivity, V = (4/3)7rr^ or , e is the 
core volume, is the permeability, ds is the character- 
istic length scale on which the field vector changes. If 
consider that 77 ~ 1 m 2 /s, r core ~ g?b and \i ~ 1, in the 
case D ~ W 1 ~ 1 TW we obtain the value of toroidal 
magnetic field Be ~ 0.3 T, which is in good agreement 
with theoretical estimations [6]. 

At the same time, in spite of the fact that the axion 
mechanism of solar dynamo-geodynamo connection ex- 
plains well the strong negative correlation between the 



magnetic field of the solar tachocline zone and the Earth 
magnetic field from the physical standpoint, it can not 
explain other correlations in Figure 1 (between the mag- 
netic field of the solar tachocline zone and variations of 
the Earth angular velocity, average global ocean level 
and the number of large earthquakes with the magni- 
tude M>7) from the energy standpoin. However, under 
certain conditions, i.e. within the framework of the hy- 
pothesis of natural nuclear georeactor existence on the 
boundary of the liquid and solid phases of the Earth 
core [16, 17], the axion mechanism can effectively pro- 
vide these correlations 



IV. SOLITON-LIKE NUCLEAR GEOREACTOR 
AND AXION MECHANISM OF THE EARTH 
CORE "HEATING" 

Now it is obvious that the magnificent experiments of 
the KamLAND-collobaration over the last 8 years [18] 
have been extremely important not only for observation 
of reactor antineutrino oscillations. They make it pos- 
sible for the first time to verify one of most vivid and 
mysterious ideas in nuclear geophysics the hypothesis of 
natural nuclear georeactor existence (see [16] and refs. 
therein). In spite of its singularity and long history, this 
hypothesis becomes especially attractive today because 
it enables to explain clearly from the physical stand- 
point different unrelated, at the first glance, geophysi- 
cal anomalous phenomena whose fundamental nature is 
beyond doubt [19]. 

We have to note that, in spite of the fact that the 
experimental KamLAND-data are well described within 
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FIG. 5. Prompt event energy spectrum of u e candidate events 
(the years 2002-2009) [17]. The shaded background and 
geoneutrino histograms are cumulative. Statistical uncertain- 
ties are shown for the data; the violet band on the blue his- 
togram indicates the event rate systematic uncertainty within 
the framework of the georeactor hypothesis. The total geore- 
actor power is 29.7 TW. Georeactors are at a distance of 6400 
and 6830 km from the KamLAND-detector. 

the framework of georeactor model [16, 17] (see Figure 5) 
and the location of soliton-like nuclear georeactors (Fig- 
ure 6 [17]) is determined by triangulation of the Kam- 
LAND [18] and Borexino data [19], some geophysicists 
have doubts not only about existence of the georeactor, 
but in the first place about its power. In this connection 
we would like to pay attention to the strange restriction 
(W < 6.2 TW) on the value of nuclear georeactor thermal 
power W, which, unfortunately, has been frequently met 
in the scientific literature recently [19-22]. This restric- 
tion terrifically masks and distorts clear understanding 
of the problem of georeactor existence, which is intricate 
enough by itself. 

Indeed, one of the conclusions of the KamLANDcollob- 
oration is the upper bound of nuclear georeactor thermal 
power (W < 6.2 TW at 90% C.L.), which is a direct con- 
sequence of uncertainty of KamLAND experimental data 
[20]. However, it is necessary to keep firmly in mind that 
this restriction is true only for the concrete parameters of 
mixing (Am^ = 7.58-10~ 5 eV 2 ,t&n6 12 = 0.56) obtained 
within the framework of the concrete x 2 -hypothesis of 
KamLAND-expcriment which takes into account the ex- 
istence of georeactor within the framework of nonzero 
hypothesis [21], but absolutely ignores such a nontriv- 
ial property of the nuclear georeactor as an uncertainty 
of georeactor antineutrino spectrum, which in the case 
of soliton-like nuclear georeactor reaches ~100% [17]. As 
shown in Ref. [17], the account of this uncertainty within 
the framework of maximum likelihood function leads (in 
the minimization of the x 2 -function) to considerable ex- 
pansion of restriction on the nuclear georeactor heat 
power (~30 TW) and, accordingly, to the new oscillation 
parameters (Am^ = 2.5 • 10~ 5 eV 2 ,tan0 12 = 0.437) for 
reactor antineutrino. 

However, in spite of obvious attractiveness of the hy- 



pothesis of natural nuclear georeactor existence there are 
some difficulties for its perception predetermined by non- 
trivial properties which georeactor must possess. At first, 
natural, i.e. unenriched, uranium or thorium must be 
used as a nuclear fuel. Secondly, traditional control rods 
are completely absent in the reactivity regulation system 
of reactor. Thirdly, in spite of the absence of control 
rods a reactor must possess the property of so-called in- 
ner safety. It means that the critical state of the reac- 
tor core must be permanently maintained in any situa- 
tion, i.e. normal operation of the reactor is automatically 
maintained not as a result of operator's activity, but by 
virtue of physical reasons-laws preventing the explosive 
development of chain reaction by natural way [24]. Fig- 
uratively speaking, the reactor with inner safety is the 
"nuclear installation which never explodes" [25]. 

It seems to be strange, but reactors satisfying such un- 
usual requirements are possible in reality. For the first 
time the idea of such a self-regulating fast reactor (so- 
called mode of breed-and-burn) was expressed in a gen- 
eral form by Russian physicists Feynberg and Kunegin 
[26] and relatively recently "reanimated" as an idea of 
the self-regulating fast reactor in traveling-wave mode of 
nuclear burning by L. Feoktistov [27] and independently 
by Teller, Ishikawa and Wood [28]. 

The discussed nuclear georeactor located on the 
boundary of the liquid and solid phases of the Earth core 
has an unique and important for our aims feature, which 
consists in the fact that the fission cross-section of 239 Pu 
(generated due to the georeactor operating) is the sharply 
nonlinear function of temperature in the range 3000-5000 
K (Figure 7a) . It means that the variations of the Earth 
core temperature generated by the mechanism of solar 
dynamo-geodynamo connection will induce correspond- 
ing variations of the nuclear georeactor thermal power. 
It's strange, but it is true and it is confirmed by inverse 
correlation between the solar magnetic field (Figure 1) 
and the nuclear georeactor thermal power over the pe- 
riod 2002-2009 (Figure 7b). Thus, on the one hand, such 
a coordinated behavior of the solar magnetic field and 
the nuclear georeactor thermal power is an indirect con- 
firmation of reality of the " axion mechanism- nuclear geo- 
reactor" energy chain and, on the other hand, according 
to the estimated power of various geophysical processes 
[29] , such a generalized mechanism can provide the solar- 
terrestrial correlations shown in Figure 1 effectively. 

If the georeactor hypothesis is true, the fluctuations of 
georeactor thermal power can influence the Earth's global 
climate in the form of anomalous temperature jumps in 
the following way. Strong fluctuations of the georeactor 
thermal power can lead to partial blocking of convection 
in the liquid core [11] and change of angular velocity of 
liquid geosphere rotation, thereby, by virtue of the con- 
servation law of Earth's angular moment to change of 
angular velocity of the mantle and the Earth's surface 
respectively. It means that the heat or, more precisely, 
dissipation energy caused by friction of earthly surface 
and bottom layer can make a considerable contribution 
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FIG. 6. Distribution of geothermal power density on the Earth [23] superposed with the conjugate "pseudogeoreactor" ellipsoidal 
closed curves, which were built on basis of KamLAND (red lines) and Borexino (blue lines) experimental data [17]. (*) operating 
nuclear georeactors; (O) an d (°) nuclear georeactors, whose power (if they are operating) is an order of magnitude or more 
less than the thermal power of reactors designated by (*). 




FIG. 7. (a) Dependence of the 239 Pu fission cross-section averaged over the neutron spectrum on fuel medium temperature for 
limiting energy (3kT) of the Fermi and Maxwell spectra. The similar dependence for the 235 U fission cross-section is shown for 
comparison, (b) Time evolution of the nuclear georeactor thermal power W. 



to the total energy balance of the atmosphere and thereby 
significantly influence on the Earth global climate evolu- 
tion [10, 11]. 



allowance for nonlinear friction by Gilmorc [30] and the 
climatic potential (7) has the form of the van der Pole- 
Duffing type equation (sec Figure 2): 



V. BIFURCATION MODEL OF THE EARTH 
GLOBAL CLIMATE ON THE ANNUAL TIME 
SCALE 

Newton's second law for friction rough surfaces (the 
Earth surface and the atmosphere surface layer) with an 



rax = — fi {x 2 — A) x — d x U, (13) 

where x is the average shift length of the atmosphere sur- 
face layer relative to the Earth surface, m is the effective 
mass of boundary layer, \i and A are parameters, U is the 
climatic potential of (7) type. 
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Using the substitutions x — Aw ■ R ■ At. £ = ujR, 
v = Auj/lj, we can write down (13) in the following form 



(At) 2 v 2 - A 



i au 



(14) 



where u is the angular velocity of the Earth rotation and 
Aw is its change over At= 1 year, R is the average Earth 
radius, v is the dimcnsionless quantity which describes 
by definition the Earth rotational velocity [7] (Figure 1). 

Since temporal variations of the global ocean level and 
temporal variations of the Earth average temperature 
strongly correlate without time lag, whereas the tem- 
poral variations of the global ocean level and temporal 
variations of the Earth rotational velocity strongly cor- 
relate with the lag ti ag ~5 years (see Figure 1), (14) with 
consideration of the approximate equality v ~ kT t -t lag 
can be rewritten in the following form 



£fcmT — —fJ-£k 



ek 2 (Atyr'-x 



T 



1 dU 
£fcAt dT ' 



(15) 



that explicitly takes into account the mechanism of solar 
power pacemaker. 

Nontrivial properties of the basic equation of bifurca- 
tion model of the Earth global climate on the annual 
time scale are exhibited in Figure 8 by the variety of 
phase portraits depending on the governing parameters 
(a, 6). Moreover, the change of shape of the assembly- 
type catastrophe potentials (7) on the plane (a, b) directly 
specifies the conditions of " warm-cold" phase transitions 
in the climatic self-oscillatory system of the van der Pole- 
Duffing type (13). 

Here it is interesting to note the following remarkable 
fact. It was found that the low order dynamic models 
of the time evolution of the toroidal magnetic field of 
the Sun derived from mean field dynamo theory are also 
described by the nonlinear oscillator equations of the van 
der Pole-Duffing type [31, 32]. In this sense, the identical 
type of equations describing the time evolution both of 
the Sun magnetic field and the Earth global climate is 
one more confirmation and, at the same time, natural 
consequence of physical, i.e., really existing, mechanism 
of solar dynamo-geodynamo connection. 

Now we return to the problem of taking into account 
of the mechanism of solar power pacemaker within the 
framework of the bifurcation model of the Earth global 
climate on different time scales. It is known that on the 
large time scales (from several to ten thousands years) on 
which our bifurcation model was considered above, the 
equilibrium state of the global climate is reached at every 



time point. It is obvious, that in this case the left-hand 
side of (15) can be set equal to zero, and (7) itself can be 
written down in the following form 



d T U' = T 3 + a'T + b' = 





a 


/ / y 


b 




FIG. 8. A plane of the parameters (a, b), the typical shapes 
of the assembly-type catastrophe potential (red lines) and the 
phase portraits (black lines on the pink squares) of the self- 
oscillating system of the van der Pole-Duffing type (13) at 
m, A = 1. Blue circles are points, to whose coordinates the 
phase portraits and the assembly-type catastrophe potential 
correspond. 



It means that the bifurcation model of the Earth global 
climate on the ten thousandth time scale really takes 
into account not only the laws of atmospheric physics, 
in particular, the laws of geometrical optics of climatic 
billiards, which generalize the cosmic rays-clouds effect 
by Sven-smark, the first (the Twomey effect) and second 
indirect aerosol effects [10], but also the mechanism of so- 
lar power pacemaker, which was masked before [10, 11] 
by renormalization procedure of the governing parame- 
ters to take into account the initial conditions. In other 
words, theoretical solutions of the bifurcation model of 
the Earth global climate on the ten thousandth time scale 
with respect to the temperature and the global ice vol- 
ume not only take into account the mechanism of so- 
lar dynamo-geodynamo connection but, in combination 
with high quality of description of the known experimen- 
tal trends of the temperature and the global ice volume 
[10, 11], are reliable confirmation of correct and holis- 
tic understanding of the basic foundations of nonlinear 
physics of the Earth global climate formation. 
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